Response of Single Polymers to Localized Step Strains 
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In this paper, the response of single three-dimensional phantom and self-avoiding polymers to 
localized step strains are studied for two cases in the absence of hydrodynamic interactions: (i) 
polymers tethered at one end with the strain created at the point of tether, and (ii) free polymers 
with the strain created in the middle of the polymer. The polymers are assumed to be in their 
equilibrium state before the step strain is created. It is shown that the strain relaxes as a power-law 
in time t as t~ v . While the strain relaxes as 1/t for the phantom polymer in both cases; the self- 
avoiding polymer relaxes its strain differently in case (i) than in case (ii): as £-( 1 +")/( 1 + 2 ") an d as 
£- 2 /( 1 +2i') reS p ec tively. Here v is the Flory exponent for the polymer, with value « 0.588 in three 
dimensions. Using the mode expansion method, exact derivations are provided for the 1/t strain 
relaxation behavior for the phantom polymer. However, since the mode expansion method for self- 
avoiding polymers is nonlinear, similar theoretical derivations for the self-avoiding polymer proves 
difficult to provide. Only simulation data are therefore presented in support of the £ - ( 1+l ')/( 1 + 2l/ ) an d 
the i _2 /( 1 + 2l/ ' behavior. The relevance of these exponents for the anomalous dynamics of polymers 
is also discussed. 
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PACS numbers: 36.20.-r, 82.35Lr, 02.70.Uu 



I. INTRODUCTION 

If a polymer is subjected to local step strain, i.e., a 
small part of a polymer is made to undergo a relatively 
fast conformational change, during subsequent evolution 
the polymer will readjust itself in an attempt to relieve 
its strain. The local conformational change will alter the 
polymer's local chain tension; and the new chain ten- 
sion will be unable to maintain the polymer in equilib- 
rium. In response to that, monomers will be pulled from 
(or pushed away to) the adjacent part of the polymer, 
thereby spreading the effect of the local strain. In time, 
the effect of the local strain will spread through the en- 
tire polymer along its backbone, before equilibrium con- 
ditions can be finally restored. 

Studies on strain relaxation in collective polymeric sys- 
tems are abundant in traditional polymer physics, such as 
for (dilute/semi-dilute) polymer solutions and for poly- 
mer melts [3]. From this perspective, how a single poly- 
mer relieves its local step strain may seem to be a purely 
theoretically motivated problem. However, experimen- 
talists' ability to manipulate polymeric systems at sin- 
gle polymer level — specially in the context of biological 
polymers, or biopolymers — have rapidly grown in the 
last few years; e.g., DNA separation in nanochannels Q, 
dynamics of RNA polymerase pj], biopolymer transloca- 
tion 0, H, 0, 0, S, [§], packagin g an d ejection of bacte- 
riophage DNA during infection jToL fllj ] su rface desorp- 
tion of polymers using a pulling force [12|. Such single 
polymer experiments have been continuously challenging 
polymer theorists; one can almost claim that polymer 
physics at a single polymer level is being reborn through 
these recent developments. Indeed, our motivation to 
study the response of single polymers to localized step 
strains, stem from the fact that there are systems whose 
dynamics are determined by the polymers' local strain re- 



laxation mechanism. Take for example polymer translo- 
cation, where the polymerpasses through a narrow pore 
in a membrane @, A translocating polymer 

is composed of two polymer strands (labeled A and B re- 
spectively), one on each side of the membrane. The only 
way the two strands interact with each other is through 
the pore: as the monomers translocate, they leave one 
strand to join the other. Monomers leaving strand A lo- 
cally increases the chain tension of strand A at the pore, 
and as they join strand B across the membrane, they re- 
duce the chain tension of strand B, also locally at the 
pore. How the segments relieve these local strains de- 
termines the dynamics of translocation fljjL [13 . HR [lR |. 
Similarly, in the case of polymer adsorption on a rigid 
surface, when a monomer gets adsorbed, it creates a lo- 
cal (at the adsorbing surface) step strain in the polymer, 
and the adsorption kinetics is governed by how the poly- 
mer relieves this strain [l7j |. 

The fact that local step strain relaxations of a poly- 
mer is governed by a power-law in time can be argued 
on general theoretical grounds. Let us consider the ap- 
plication of the step strain of magnitude eo at a given 
location (say the n*-th monomer) of a polymer of length 
N at t = 0. This strain will excite all fluctuation modes 
of the polymer. The amplitude a q of the q-th mode 



tpq can be obtained from the equation eo = 



,(°) 



q = 1,2, ...N, Typically, in polymer physics, the q-th 
fluctuation mode of a polymer has an associated relax- 
ation time T q ~ (N/q)P for some (3, where rjy ~ is 
the longest relaxation time of the polymer, correspond- 
ing to the slowest mode q = 1 of the polymer ((3 = l + 2v 
for a Rouse polymer, and j3 = 2>v for a Zimm polymer). 
The subsequent evolution of this strain will then be given 

by e(t) = J2 q a g°Vg ex p(— V T g)- The local contribution 
of these summed over large number of exponentials at 
n* will yield a power-law, implying that e n * (t) ~ t~ v for 



2 



some T), multiplied by the overall terminal exponential 
decay ~ exp(— i/rjv)- Such power-laws are often referred 
to as "memory effects". The quantity 77, the exponent for 
the power-law, characterizes the response of single poly- 
mers to local step strains. For the two physical systems 
discussed above, namely polymer translocation and ad- 
sorption of polymers on rigid surfaces, it is the exponent 
77 that dictates the dynamics fljjL \\A. KEl 

The purpose of this paper is to report the exponent 77 
for phantom and self-avoiding polymers in three dimen- 
sions in the absence of hydrodynamic interactions. The 
specific way we create the local strain in the polymers is 
as follows. At a given location (say the n*-th monomer) 
of an equilibrated polymer of length N, we inject p (<C N) 
crumpled monomers at t = 0, bringing its length to N+p. 
Following the monomer injection at t = 0, apart from the 
newly injected monomers, the polymer follows random 
walk (or self-avoiding walk) statistics, i.e., the strain in 
the polymer is localized at monomer number n*. In the 
subsequent evolution of the polymer, we then keep track 
of how these p crumpled monomers unfold themselves, 
which yields us the exponent 77. Note that the specific 
way we choose to create the local strain in the polymers 
is indeed motivated by the actual microscopic dynam- 
ics of polymer translocation or polymer adsorption on a 
rigid surface: as remarked above, for polymer translo- 
cation it is the addition or disappearance of monomers 
to the polymer segments on either side of the membrane 
that creates the local strain (and similarly for the case of 
polymer adsorption on a rigid surface). 

We calculate 77 for two different cases each for three- 
dimensional phantom and self-avoiding polymers: (i) 
polymers tethered at one end with the strain created at 
the point of tether, and (ii) free polymers with the strain 
created in the middle of the polymer. We derive that 
n = 1 in both cases; however, for the self-avoiding poly- 
mer we show that 77 = (1 + v)/(l + 2u) for case (i), and 
77 = 2/(1 + 2v) for case (ii). Here v is the Flory exponent 
for the polymer, with value ~ 0.588 in three dimensions. 
We provide exact derivations for the 1/t strain relaxation 
behavior for the phantom polymer using the mode ex- 
pansion method. The mode expansion method for a self- 
avoiding polymer is nonlinear, and hence similar theoret- 
ical derivation for 77 for the self-avoiding polymer proves 
difficult to provide. Only high-precision simulation data 
are therefore presented in support of the £ _ ( 1 + I/ )/( 1 + 2l/ ) 
and the i _2 /( 1 + 2l/ ) step strain-relaxation behaviors of the 
self-avoiding polymer. 

Although the problem of local step strain relaxation 
behavior in the polymers is motivated in this paper in 
view of polymer translocation and polymer adsorption, 
note that both physical processes correspond to the case 
(i) while the tether point lies on a rigid surface. The 
presence of the surface, in principle, can influence the 
strain relaxation mechanism, and alter the value of 77 from 
its value in the absence of the surface. However, since in 
Refs. [H, 0, [HI it was shown — using a model that 
allowed direct observation of the local strain relaxation 



— that 77 = (1+v) / (l+2v) for a self-avoiding polymer for 
the case of (i) in the presence of a rigid surface as well, the 
result of this paper therefore implies that the local strain 
release mechanism for self-avoiding tethered polymers is 
unaffected by the presence of a surface at the tether point. 
Note that recently, albeit indirectly, a different polymer 
model has confirmed that 77 = (l + v)/(l + 2v) for a self- 
avoiding polymer for the case of (i) in the presence of a 
rigid surface [H, [H|, in support of Refs. [ll [II [13 . 

This paper is organized as follows. In Sec. Ill Al we use 
the mode expansion technique for a phantom polymer for 
the case of (i) and derive that 77 = 1. In Sec. Ill Bt we 
then consider case (ii) for a phantom polymer to again 
derive that 77 = 1. In Sec. IHII we report the corresponding 
results for self-avoiding polymers, and discuss the reasons 
why the self-avoiding behaves differently in case (i) than 
in case (ii). The paper is then concluded in Sec. IIVI with 
a discussion on the relevance of these exponents for the 
anomalous dynamics of polymers. 



II. RESPONSE OF PHANTOM POLYMERS TO 
LOCAL STEP STRAIN 

With r(n,t) as the physical location of the 77,-th 
monomer of the polymer at time t, we start with the 
Rouse equation for a phantom polymer and add thermal 
noise f(n, t) to it: 



dr 
di 



(1) 



In Eq. (HJ) the thermal noise f(n, t) satisfies the property 
that (f(n,t)) = and {/ a (n,t)/ /3 (n / ,t')> = 25 a/3 S(n - 
n')S(t — t'); a, (3 = x, y, z. For case (i), the polymer with 
its zeroth monomer tethered at the origin we define the 
g-th mode for a polymer of length (N+p), tethered to a 
fixed point at the origin as [l| 



x g (t) = 



N+p Jo 



N+p 



dn svn.(k q n) f(n, t), 



(2) 



with k q = 2(]y~+ y an< ^ ^ = 1 1 2, 3, ... , and similarly f q , 

the q-th mode for the thermal noise. The sine-expansion 
in Eq. ^ satisfies the boundary condition that f(0, t) = 

dr(n, t) 



Vt, and also that at the free end 



dn 



= 0. For 



N 



case (ii) we define the g-th mode for a polymer of length 
(N + p), moving freely in space as [H] 



x q (t) 



1 



N+p 



with kn = , x . 

q (N+p) 



N+p Jo 
Tvq 



dn cos(k q n) f(n, t), (3) 
, and q = 0, 1, 2, 3, . . ., and similarly 



f q , the q-th mode for the thermal noise. In this case 
the cosine-expansion satisfies the boundary condition 



3 



that that at the free ends of the polymer 



df(n, t) 



df(n, t) 



dn 



for the end-tethered and free polymers respectively. 



di 



N 



In terms of the transforms ([2j) and §2§ the Rouse equa- 
tion {!]) reduces to the Langevin form 



dX q 

at 



-k 2 q X q + f, 



q • 



(4) 



where f q is defined similar to Eq. ([2]) [resp. Eq. (|3jl]. 
This reduction to the Langevin form also yields 



(/pa(*)) = 5 </ P aW/^(*')) = 



1 



N+p 



3pq $a(3 



In terms of X q (t) the monomer locations in physical space 
are then given by 



? (n, t) = 2 s'm(k q n) X q (t) and 
i 

i,t) = 2^2 cos{k q n) X q (t). 



r n, 



(6) 



A. Local strain relaxation for case (i): 
end-tethered phantom polymers 



As we crumple the extra p (<C N) monomers at the 
tether point to an equilibrated polymer of length N at 
time t = 0, the length of the polymer instantaneously be- 
comes N + p. The ensuing time-evolution of the polymer 
is then described by 



X q (t) = e- k « t X q (0) + / dt' e-M*-*') f q { t >) , ( 7 ) 
Jo 



i.e., 



f(n, t) = 2 \ sm(k q n) 



1 



-hi t 



X q (0) + / dt' (*-*') f q (t>) 



(8) 



After the injection of p monomers at t — 0, to follow the deviation from random-walk statistics along the polymer's 
backbone at a given location of the polymer, say at monomer number no, we consider another nearby monomer m, 
define n = \n\ — no\ and r 2 (n, t) = [f(ni, t) — r(no, t)] ■ [f(ni, t) — r(no, t)] 



r 2 (n,t) = ^ [sin(fc g ni) - sin(fc 9 n )] 

A 9 («l,Tio) 



[sin(fco'^i) — sin(fc 9 'no)] 



e- fe ?*J? g (0)+ / dt' e'^^-^ f q {t') 
Jo 

- k l> 1 X q , (0) + / dt" e~ k l' ( *- t " ) £ (t") 



(9) 



A i (ni,n ) 



With the aid of Eq. J5]), for a given polymer realization at t = 0, the average over the evolution histories (i.e., noise 
realizations) for t > 0, denoted by the angular brackets (.), for this polymer yields 



(r 2 (n,t)) =4^ ^A q (n 1 ,n )A ql (n 1 ,n ) , 



-(K+k 2 ,)t 



X, 



q.q' 



,(0)-AV(0)]} 



6 



(N+p) 



A2 q (. n li n o) 



i - e - 2fe ;* 



(10) 



At t — > oo, the t-dependent terms drop out, leaving us with 

6 \ [sin(fc g ni) — sin(fc 9 no)] 2 6 



(r (n, t — ► oo)) 



(N + p) 



1-2 
'/ 



[sin(niz) - sin(n a;)] 2 _ 

Ct3j ~ — tjiL . 



(11) 



which confirms that the polymer returns to equilibrium as t — > oo, as it should. 

Since the strain at t = is created at the tether point, i.e., at monomer number zero of the polymer (of length 
N+p), to quantify its relaxation we track (||r 2 (n,t)||) by choosing n Q = n* = and n x — n, with n ~ 0(p). Here ||.|| 
denotes a second average over equilibrated configurations of the polymers at t = 0. From Eqs. ([9]) and (JH), we can 



4 



then write 

(||r 2 (M)||H3n + 4i:{sin(fc^^ (12) 



<]■'( 



Notice that if the polymer of length (N + p) were already at equilibrium at t = (i.e., no step-strain were created 
anywhere in the polymer), then it would have remained in equilibrium Vi > 0; i.e., (||r 2 (n,i)||) = (||r 2 (n,i)||( eq )) = 
3n Vi. In that case, Eq. ffl2|) would reduce to 

4]T {sin(fc g n) sm(Vn) ||X^>(0) ■ X^(0)||} = — ° £ e"^*, (13) 

9,9' 9 9 

where X q eq '(0) is obtained from Eq. ([2]) for the polymer at equilibrium at t = 0. An explicit calculation of Eq. lfT3|) 
has also been provided in Appendix A [Eqs. (|AHIA5[) ] . 

Based on Eq. (fl3|) we can now replace the last term on the r.h.s. of Eq. lfT2l) by the l.h.s. of Eq. (fl3l) to write 

(| \r 2 {n, t) 1 1) - 3n = 4 ^sin(/c 3 n) sin(fc^n) e'^+'v >* <? g , g ,, (14) 

9,9' 

with g q . q . = \ \X q (0) ■ X ql (Q)\\ - | |^ eq) (0) • A<f q) (0)||. The quantity gfl has already been simplified in Eq. ^ as 



„< 2 > 



f,,,' 9 



(2) _ 3 1 x n 

9 ™' _ (N+p)2k q k q . [ j 

while the quantity is explicitly evaluated in Appendix B [Eqs. (|B1IIB4[) ] . Having combined these two quantities, 
in the limit of p — > we find that 

ra "(iv + 'kv (16) 



which, when used in conjunction with Eqs. lfT2[) and (|T4|l . we obtain 

-i 2 



|r (n,tj||)-3n + ^ ^2 

9,9' 4 4 



saxinx) e 
ax 



x 



3n-— (17) 



at long times. In other words, the local strain at the tether point relaxes as 1/i; i.e., the local step strain relaxation 
exponent r) = 1. 

B. Local strain relaxation for case (ii): free phantom polymers 

For the local strain relaxation following the injection p crumpled monomers at n* = N/2 into freely moving phantom 
polymer at t — we follow the same route as in Sec. Ill A} however, one needs to replace the sine-expansion by cosine- 
expansion. While Eqs. |7l fT2|) are trivially reproduced with this replacement, for the rest of the calculation we need 
two small modifications. The first one of them is to choose ri\ = (N + p — n)/2 and n = (N + p + n)/2 such that 
(||r 2 (n,t)||), as defined above Eq. ((9|), can once again quantify the local strain relaxation of the polymer. The second 
one is that A q (nx,no) is now defined as A g (ni,no) = [cos(fe g ni) — cos(fe g no)]. These lead us to the equivalent forms 
of Eqs. (gHH as 



r 2 (M)IIH 3n + 4^(n 1; no)^!, no) { 



(N+p) 

9,9' 9 



and [as explicitly evaluated in Eqs. l|A6llAlip in Appendix A] 

9,9' '9 q 



5 



Similarly, analogous to Eq. (fT4| we have 

(| |r 2 (n, t) 1 1) - in = 4 ^ A 9 (m , n )^ (m , no) e - ( *W' g q . q , , (20) 

where = ||*,(0) • X q ,(Q)\\ - ||X^(0) • xj q) (0)||, with 
-> ' >. / 



(2) = O J_ r 

(iV+p)2fc 9 fc g , fe "' fe «' 



(21) 



from Eq. (fl9|) . The explicit evaluation of g q L q , is carried out in Appendix B [Eqs. (|B5HB7|) ], Having combined g^, 

(2) 

and g ',, below we present the final result for g qjq > in the limit of p — > 0: 



9q,q' 



3p 



iin[k q (N + p)/2] sm[k q , (N + p)/2] 



(iV+p) 2 

which, when used in conjunction with Eq. lf20|) . we obtain 

<lk 2 (M)IIH3n- 



(22) 



12p A g (ni, n )A g >(ni, Ttp) sin[fc g (7V + p)/2] sin[fc g /(iV + p)/2] e 

2 / /-i- ^ A > 



(N+p) 



kq kq' 



Finally, with A q {nx, no) = 2sin[fc 9 (iV + p)/2] sin[fc 8 n/2], and sin[fc g (iV +p)/2] = sin[7rg/2] for q = 1,2,3,..., Eq. JH 
reduces to 

<||r 2 (M)||)=3n 



(iV + p) 2 



, sin(fc g n) sm 2 [fc g (iV + p)/2] e" fe '* 



-1 2 



L q 



3n — 



(N+p)< 



E 

q € odd 



sin(k q n) e « 



= 3n 



24p 



rf.7' 



sin(nx) e 



3n-*£ (24) 



which, just like Eq. ifTTj) . approaches its asymptotic value 3n as 1/t; i.e., once again the local step strain relaxation 
exponent rj = 1. 



III. RESPONSE OF SELF-AVOIDING 
POLYMERS TO LOCAL STEP STRAIN 



We use a Monte Carlo based lattice polymer model 
to study the local step-strain relaxation for self-avoiding 
polymers. In this model, the polymer consists of a se- 
quential chain of monomers, living on a FCC lattice. 
Monomers adjacent in the string are located either in 
the same, or in neighboring lattice sites. Multiple occu- 
pation of lattice sites is not permitted, except for a set 
of adjacent monomers. The polymer moves through a 
sequence of random single-monomer hops to neighboring 
lattice sites. These hops can be along the contour of the 
polymer, thus explicitly providing reptation dynamics. 
They can also change the contour "sideways", providing 
Rouse dynamics. The reptation as well as the sideways 
moves are attempted with rate unity, which provides us 
with a definition of time in this model. This model has 
been used before to simulate the diffusion and exchange 



of polymers in an equilibrated layer of adsorbed poly- 
mers [20||, polymer translocation under a variety of cir- 
cumstances (la . [ill. [2l[ , and polymer adsorption 
to rigid surfaces [IT]]. Multiple occupation of the same 
site by adjacent monomers of the polymer, in this model, 
gives rise to "stored lengths" (see Fig. 2 of Ref. for 
an illustration). Upon injection of p extra monomers into 
the polymer at the lattice site where the n*-th monomer 
[n* = and N/2 for cases (i) and (ii) respectively] is 
located at t = the local stored length density is im- 
mediately increased by p. To measure the local strain 
relaxation of the polymer we therefore track the density 
of stored lengths per monomer in these new p monomers, 
p p {t) as a function of time. Of course p p (t) would ap- 
proach some "offset" value po as t — > oo. 

We have already argued in the introduction that the 
strain-relaxation behaves as t~ v exp(— t/r N ). The ter- 
minal exponential decay exp(t/r N ) with t n ~ N 1+2v 
is expected from the Rouse relaxation dynamics of the 
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entire polymer. To understand the physics behind the 
exponent 77, we use the well-established result for the 
relaxation time t n for n self-avoiding Rouse monomers 
scaling as t n ~ n 1+2v . On the basis of the expres- 
sion of t n , we anticipate that following the injection of 
p monomers at t = 0, by time t the extra monomers will 
be well-equilibrated across the inner part of the poly- 
mer up to n t ~ t 1 /( 1+2,y ) monomers around n* , but not 
significantly further. This internally equilibrated sec- 
tion of (n t + p) monomers extends only to r(n t ) ~ n v t , 
less than its equilibrated value (nt + p) v , because the 
larger scale conformation has yet to adjust to the lo- 
cal strain. As a result, internally equilibrated section 
of (nt +p) monomers remains at a state of excess free en- 
ergy 8F ~ kBT[6r(nt)/r(nt)] 2 . The excess p monomers 
need to find their own physical space by pushing the other 
monomers away for both cases (i) and (ii) , but for case (i) 
as the zeroth monomer remains tethered, we expect them 
to feel a force of magnitude / derived from the excess 
free energy as / = dF/dr(n t ) ~ fcgT 5r(nt)/r 2 (n t ) ~ 
£-(i+«0/(i+2i') j wn i cn dictates the relaxation of the step- 
strain; i.e., f] = (1 + v)/(l + 2v). In case (ii) however, 
the force derived from the excess free energy does not 
yield 77, as the internally equilibrated section will simply 
move under the effect of the force. Instead, in case (ii) 
we expect these p monomers to feel a chemical poten- 
tial of magnitude /i derived from the excess free energy 
as p, = dFjdn t = [dF/dr{n t )][dr(n t )/dn t ] ~ £-2/(1+2*). 
The step strain relaxation is then dictated by the chemi- 
cal potential p.; i.e., n = 2/ (l+2v). In Fig. [IJ by tracking 
p$(t) for N = 195 and p — 5, we provide confirmation 
of this physics. Note that the result for r\ for case (i) is 
consistent with the corresponding two dimensional case 
in Ref. [Hjj |. as it should be. 



IV. DISCUSSION 

In this paper, response of single polymers to localized 
step strains is studied for two cases in the absence of 
hydrodynamic interactions: (i) polymers tethered at one 
end with the strain created at the point of tether, and 
(ii) free polymers with the strain created in the middle 
of the polymer. The polymers are assumed to be in their 
equilibrium state before the step strain is created. Using 
mode expansion technique for Rouse equation it is shown 
that for phantom polymers in both cases the strain re- 
laxes in time as 1/t. However, for self-avoiding polymers 
for the two cases the strain relaxes as £~( 1 +< y )/( 1 + 2,y ) and 
as £~ 2 /( 1+2l/ ) respectively. The strain relaxation behav- 
ior £ _ ( 1 + i/ )/( 1 + 2l/ ) for a self-avoiding polymer for case (i) 
is consistent with an earlier reported result in two di- 
mensions [15J. Based on the results reported here, and 
combined with those of Refs. [T3 . we can conclude 

that the result for case (i) is independent of the presence 
of a surface at the tether point. 

Although in both cases (i) and (ii) the local step strain 
puts the polymer in a state of excess free energy, the dif- 



10° 




t 



FIG. 1: (color online) Numerically differentiated data for 
\dp 5 (t)/dt\ for cases (i) [top set of points, in red] and (ii) [bot- 
tom set of points, in blue], for N = 200 and p = 5 (10, 000, 000 
realizations each), showing the respective t~^ 1+v ''^ +2v ' > (top 
straight line, in red) and t~ 2 ^ 1+2v ^ (bottom straight line, in 
blue) power-law decay for ps(t). Note that (1 + v) / (1 + 2u) w 
0.73 and 2/(1 + 2u) w 0.92. We use numerical differentiation 
in order to remove the t — > 00 offsets of ps(t). The data for 
case (i) is displaced upwards by a factor 2 in the y-direction. 
Inset: Ratio r(t) of the \dp*,(t)/dt\ values for cases (i) and (ii), 
showing that r(t) follows the power-law t( 1 ~">/( 1 + 2l/ ) ■ where 
the value of (1 — u)/(l + 2u), the difference in the values of rj 
for cases (i) and (ii), is « 0.19. 

ference between the results for the self-avoiding polymers 
for these two cases stems from the fact that the tether 
point provides a point of reference for the polymer in 
case (i), but not in case (ii). As a result, for case (i) we 
need to consider the force, while for case (ii) we need to 
consider the chemical potential, derived from the excess 
free energy. For phantom polymers however, since differ- 
ent parts of the polymer do not interact with each other, 
there is no need for the strained monomers to physically 
push away the other monomers of the polymer in order 
to be able to relieve their strain, and hence for case (i), 
the force derived from the excess free energy plays no 
role in the localized strain relaxation for the phantom 
polymer. In fact, precisely because of the same reason, 
we expect to see 1/t strain relaxation for phantom poly- 
mers also in the presence of a surface at the tether point. 
With t^ 1 = i^ 2 /( 1+2,y ) for phantom polymers (v = 0.5), 
the relevance of this paper is that one cannot trivially 
extend the local strain relaxation behavior for tethered 
phantom polymers to self-avoiding polymers by replacing 
v = 0.5 by v ~ 0.588 in three dimensions. 

In earlier published works [3, 0, EE EE H3], a 

"voltage-current" relationship <p(t) = \ dt' p(t — t')s(t') 

Jo 

between s(t), the instantaneous rate of translocation, and 
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the polymer's chain tension imbalance <p(t) across the 
pore was established, where fx{t) is the memory effect de- 
rived from the polymer's local strain (alternatively, the 
chain tension) relaxation behavior at the pore. Here s{t) 
is the number of the monomer located in the pore at 
time t. Using ji(t ) ~ t~( 1 + !y )/( 1 + 2,y ) for unbiased polymer 
translocation [3, EE H3| as in case (i) for self-avoiding 
polymers in this paper, the anomalous dynamics, charac- 
terized by (As 2 (t)), where As(t) is the total number of 
monomers translocated through the pore in time t, was 
then derived by using the fluctuation-dissipation theo- 
rem, where the angular brackets denote an ensemble av- 
erage. It was found that for a translocating polymer of 
length AT, (As 2 (t)) ~ t(*+ 1 0/( 1 + 2l up to tne Rouse time 
tn ~ N 1+2u , and since no memory can survive in the 
polymer beyond the Rouse time, (As 2 (t)) ~ t for t > tjv, 
i.e., the pore-blockade time scaling as N 2+u . This result 
for the scaling of the pore-blockade time is in good nu- 
merical agreement with that of Refs. [IE obtained 
using completely different polymer models. Furthermore, 
having exploited the same "current-voltage" relationship 
between s(t) and the chain tension difference <f>{t) across 
the pore and that /i(t) - H 1 +* / >/( 1 + 2 * / > for field-driven 
translocation as well, the exponent N^ 1+2u ^ '' scaling 
was later found for the pore-blockade time for field-driven 
translocation of a polymer of length N (this result 
has recently been confirmed pl| using another different 
polymer model). Similarly, for the non-equilibrium dy- 
namics of single polymer adsorption to solid surfaces, the 
adsorption time for a polymer of length N at weak ad- 
sorption energies was also found to scale as N^ 1+2v ^^ 1+v ^ 
[171 ] . These results, put together with the discussions 
in the above paragraph [namely that the value of n for 
case (i) is independent of the presence of a surface at the 
tether point], lead us to expect that the pore-blockade 
time for unbiased translocation should scale as N 2+u for 
self-avoiding polymers, and as N 2 for phantom ones, irre- 
spective of whether translocation proceeds through a nar- 
row pore in a membrane or whether it proceeds through 
a narrow ring (i.e., a pore without a membrane). 

It is imperative to ask, based on the local strain re- 
laxation result for case (ii), whether it would be possible 
to derive an expression for the mean-square-displacement 
(Ar 2 (n, t)) of the n-th monomer in physical space in time 
t, by tracking the physical location f(n,t) for the n-th 
monomer of the polymer at time t. In order to answer 



this question, let us reconsider the "voltage-current" rela- 
tionship between the chain tension imbalance across the 
pore and s(t), and note that for translocation s(t) is a 
scalar variable, while r(n,t) is a vector, and as a result, 
deriving (Ar 2 (n,i)} in a similar manner is more compli- 
cated. To illustrate this difficulty, let us return to the 
deterministic part of Eq. ([J): by first expressing r as a 
function of the polymer's contour /, and then expressing 
the I as a function of n, Eq. ([I]) reads 



df(n, t) 



d 2 f ( dl \ df d 2 l 

dt ~dP [dnj + ~dldri 2 ' 
The first term on the r.h.s. of Eq. lj25j) is a force that acts 
on the n-th monomer perpendicular to the contour of the 
polymer at the location of the n-th monomer at time t, 
while the second term is a force on the n-th monomer 
that acts along the contour. Note also that the term 

d 2 l dl 
-—^ is precisely the imbalance in the chain tension — — 
onr ' on 

at the n-th monomer. In the case of translocation, the 

fact that the motion of the monomer perpendicular to 
the polymer's contour in the pore is completely blocked 
means that the motion of the monomer in the pore is de- 
termined entirely by the chain tension imbalance across 
the pore. For a free polymer however, the first term on 
the r.h.s. of Eq. f25|) does contribute to the motion of 
the n-th monomer, but what is its precise contribution to 
(Ar 2 (n,t)) is not entirely clear. Nevertheless, if we con- 
sider the second term alone, then it does allow us to write 
a voltage-current relationship (exactly the same as that 
of Refs . [H, Hi, EE EE HI ) between the chain tension im- 
balance at the n-th monomer and the along-the-contour 
velocity component of the n-th monomer, but this time, 
following the polymer's local strain relaxation behavior 
for case (ii), with fx{t) ~ t^ 2 ^ 1+2u \ The application 
of the fluctuation-dissipation theorem would then imply 
that (Ar 2 (n,i)) should increase as t 2 /( 1+2iy ) along the 
polymer's contour, i.e., in physical space (Ar 2 (n,t)) ~ 
t 2v/(i+2v)^ tin the Rouge time Tn _ N i+2v. this is a well- 
known result in polymer physics. 
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APPENDIX A: DERIVATION OF ||X ? (eq) (0) ■ {0)\\ FOR PHANTOM POLYMERS 

Here we provide a derivation of Eq. lfT3|) for case (i) and an analogous form of it for case (ii) . 
For case (i), by definition 

1 f N+p f N+p 

||Af <%) • 4^(0)11 = j^—^ J o dn sm(k q n) J dn> S m(k q m') ||r» • f(n')l l (eq) • (Al) 



8 



In equilibrium the polymer satisfies random walk statistics along its entire backbone. Hence, with Q(x) denoting the 
Heavyside function of x, 



\\f(n) ■ r(n')|| (eq) = 3nO(n' - n) + 3n'Q(n - ri), 



(A2) 



which reduces Eq. (|A1[) to 



||XM(0)-X^(0)|| 
3 



(N+p) 2 
3 

(N+p) 2 
3 



N+p rN+p rN+p rN+p 

dnnsm(k q n) / dn' sm(k q rn') + / dn' n! sm(k q in') / dnsm(k q n) 

(i in JO Jri 

N+p sin(fc g n) cos(k q >n) r N +P / r sm(k q in') cos(k q n') 



(N+pf 



sin[2k q (N +p)} - 2k q (N+p) cos[2k q {N + p)] 



4fc3 



,K- +(l-4-„fc,) x 



kg cos[k q , (N+p)} sm[k q (N+p)} - cos[k q (N+p)} { (k 2 - k' 2 )(N+p) cos[k q , (N+p)} + k q , sm[k q > (N+p)} } 

kqkq> (k 2 — k 2 ,) 



(A3) 



The second step of Eq. (|A3[) requires cos[k q (N+p)} = cos{k q > (N+p)} = 0, while in the last step using cos[k q (N+p)} = 
cos[k q ,(N + p)} = 0, we first see that |A^ eq) (0) • A^ q) (0)|| cx 8 kq ,h ql , and moreover, with sm\2k q (N + p)} = and 
cos[2fc (? (A r + p)} = — 1, we obtain 



}x^\o)-x^m\ 



(N+p) 2k 2 



fa, 



2 : 



(A4) 



i.e., 



9,9' 



{sin(fc 9 n) sm(k q ,n) e~ {k > k > | }X^ (0) • X q ^ ] (0)| | } 



6 



(N+p) 



E 



sin (fc g n) 2fc 2 t 

1,2 
9 



(A5) 



To derive a similar expression for||x( eq) (0)-Xy q) (0)||for case (ii) we express r(n, 0), the physical location of the n-th 
monomer at t = 0, relative to r(0, 0), the physical location of the first monomer at t ~ as 0) = r(0, 0) + f'(n, 0). 
Then 



*,(0) 



1 



N+p 



dn cos(k q n) r(n, 0) 



1 



(N+p) J 



N+p 



dn cos(k q n) [f(0, 0) + r'(n, 0)], 



(A6) 



implying that 



||J?M( ).X^(Q)| 



1 



7V+p 



l r N +p 



N+p 



dn cos(k q n) / dn' cos(Ayn') [||r 2 (0,0)|| + ||f'(n) • r '(tz')||] 



(eq) 



, . dn cos(fc„n) / c?n' cos(Avn') llr '(n) • r'(n')|K eq ). 
(N + p) 2 J a Jo 



(A7) 



To obtain the second step of Eq. I|A7[) ||r 2 (0,0)|| = has been used by a trivial translation of origin to obtain 
r(0, 0) = 0, without affecting any part of the calculation. 

In terms of r'(n,0), we can once again use 



\f(n) ■ r(n')}} p= a — 3nO(n' — n) + 3n'Q(n — n'), 



(A8) 



which reduces the expression for ||x( eq) (0)-Xj q) (0)|| to 
||X( eq )(0).X^(0)|| 

N+p rN+p pN+p pN+p 

dnncos(k q n) / dn' cos(k q in') + / dn' n' cos(fc ? <n') / dncos(k q n) 

<< Jn JO Jri 

+p cos(k q n) sm.(k„in) r N +p .cos(kqin')sw.(k q ri') 
dnn : h / an n 



(N+p) 2 
3 



(N + p) 2 
3 

(N + P y 



o *V Jo 

2fc g (AT + p) cos[2k q (N + p)] - sin[2fc g (iV + p)] 
4fc3 



4„fe , +(1-Su k „,) x 



fc g > cos[fc^(7V+p)]sin[fc g (7V+p)] - sin[£y (A^+p)] {(fc^-fcf )(7V+p) sin[fc 9 (7V+p)]+fc g cos[fc g (iV+p)]} 



(A9) 



The second step of Eq. (|A9[) requires sin[fc g (A^+p)] = sin[k q '(N+p)] = 0, while in the last step using s'm[k q (N +p)} 



sin[k q '(N + p)] = 0, we first see that ||A^ eq) (0) ■ A^f q) (0)|| oc 8k q ,k and moreover, with sm[2k q (N + p)] = and 
cos[2k q (N +p)] = 1, we obtain 



ll^ eq) (0)-A^ q) (0)|| 



A 

(iV + p) 2k 2 kq ' k "" 
Equation ifT0|) then yields us 

4^^(n 1 ,no)A g ,(n 1 ,no) { e' {k > k l' ]t \\X^(Q) ■ A^ q) (0)||} 



6 A q ( ni ,no) c _ 2fc ; f 



</■</' 



(7V + p) 



£.2 



(A10) 



(All) 



APPENDIX B: DERIVATION OF g {1) , FOR PHANTOM POLYMERS 

a a. a' 



To evaluate g^, for case (i) we note that X q (0) = 



by construction, X a (0) = — — , 

N+pJ 



N 



N + p 



N+p 



dn sin(fc g n) r(n, 0), and since r(n, 0) = for n < p 



dn sm[k q (n + p)] r(n + p, 0), and hence 



9, 



(i) 



1 



AT r N 



dn / dn' sin[fc 9 (n + p)] sm[k q i(n' + p)] \\r(n + p) ■ f(n' + p)\\. 



(Bl) 



* q ' q ' (N + p) 2 Jo "''Jo 

Since the polymer was in equilibrium before the p crumpled monomers were injected at the tether point, we can write 

\\r(n + p) ■ r(ri + p)\\ = 3nQ(n' - n) + 3n'6(n - n'). (B2) 

Thereafter, using Eq. (|B~2jl . and cos[k q (N + p)] = cos[k q > (N + p)} = sin[(fc g - k q ,)(N + p)] = sin[(fc 9 + k q ,)(N + p)] = 
sm+2k q (N +p)] = Q and cos[2fc g (iV + p)] = —1, the expression for g^, in Eq. <|B1|) simplifies as 



0, 



(A^+p) 2 
3N 



JY 



JY 



JY 



dnn sin[fc g (n + p)] / dn' sin[fc g / (nf + p)] + / <in' n'sin[fe g '(n' +p)] / cin sin[fe g (n +p)] 



iY 



3 sin[(fc g + fc g >)p] 3 sin[(fc g - k q >)p] . 



2(N+p) 2 k 2 w (N + p) 2 2k q k q ,(k q + k q> ) (N +p) 2 2k q k ql (k q -k q ,) 



(B3) 



In the limit p — > the two terms proportional to 5k q ,k q , in Eq. 1|B3|1 cancel each other, as Eqs. lfT5|) and (|B3[1 then 
leave us with 
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To evaluate g^ q , for case (ii) we express r(n, 0), the physical location of the n-th monomer at t = 0, relative to 
r(0, 0), the physical location of the first monomer at t — as f(n, 0) = r(0, 0) + r '(n, 0), to obtain 



1 



9, 



(i) = 

(n + p ) 2 J 
1 



rf/? 



dn' cos(k q n) cos(k q >n')\\r'(n,Q)-r'(n',Q)\\ 



(N+p) 2 



N+p rN+p 



N+p 



dn / dn' cos(k q n) cos(k q in')f(n) + / dn' / dn cos(k q n) cos(kq'n')f(n r ) 



N+p 



(B5) 



where f(n) = [3nQ(N/2-n)+3N/2@(n-N/2)e(N/2+p-n)+3(n-p)0(n-N/2-p)].Thereaher, with sm[k q (N+p)} 
sm{k q ,(N+p)} = sm[(k q - k q ,)(N+p)} = sm[(k q + k g >)(N + p)] = and cos[2fc,(7V + p)] = 1, we find 



(i) _ 3 2k q N + 2 sm(k q p) cos[k q (N+p)] 

9„.„> - /Ar , -N9 — °ko ,k„, 



(N+p) 2 4fe3 «' (N+p) 2 

cos[(k q - k q ,)(N +p)/2] sin[(fc 9 - k q ,)p/2] cos[(fc (J + fc,,)(iV + j»)/2] sin[(fc g + k q ,)p/2] 



k q k q > (k q k q f ) 
In the limit p <C N eq. (|B6|1 can be expanded to obtain 



(i) _ 3 _ 3p 

~ 2(7V+p) fc ? fc 9 , (iV+p) 2 



kqkq' {kq ~f~ feg') 



jin[fc g (7V + p)/2] sin[fc g >(JV +p)/2] 



(1-4„ V ). (B6) 



(B7) 
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